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Quantum phase fluctuations play a crucial role in low dimensional systems. In particular they 
prevent true long range phase order from forming in one dimensional condensates, even at zero 
temperature. Nevertheless, by dynamically splitting the condensate into two parallel decoupled 
tubes, a macroscopic relative phase, can be imposed on the system. This kind of setup has been 
proposed as a matter wave interferometer, which relies on the interference between the displaced 
condensates as a measure of the relative phase between themi-. Here we show how the quantum 
phase fluctuations, which are so effective in equilibrium, act to destroy the macroscopic relative 
phase that was imposed as a non equilibrium initial condition of the interferometer. We show that 
the phase coherence between the two condensates decays exponentially with a dephasing time that 
depends on intrinsic parameters: the dimensionless interaction strength, sound velocity and density. 
Interestingly, at low temperatures the dephasing time is almost independent of temperature. At 
temperatures higher than a crossover scale T* dephasing gains signiflcant temperature dependance. 
In contrast to the usual phase diffusion in three dimensional trapa-, which is essentially an effect of 
confinement, the dephasing due to fluctuations in one dimension is a bulk effect that survives the 
thermodynamic limit. 



I. INTRODUCTION 

The existence of a macroscopic phase facilitates obser- 
vation of interference phenomena in Bose-Einstein con- 
densates (BEC). For example, an interference pattern in- 
volving a macroscopic number of particles arises when a 
pair of condensates is let to expand freely until the two 
clouds spatially overlap. Since the pioneering experiment 
by Andrews et afe, which demonstrated this effect, it has 
been a long standing goal to construct matter wave inter- 
ferometers based on ultra cold atomic gases. Such devices 
have promising applications in precision measurement^, 
and quantum information processing^, as well as funda- 
mental study of correlated quantum matter. The key re- 
quirement for interferometric measurement is determin- 
istic control over the phase difference between the two 
condensates. That is, the relative phase must be well 
defined and evolve with time under the sole influence of 
the external forces which are the subject of measurement. 
This requirement was not met for example in the setup 
of Ref. 0, where the two condensates were essentially 
independent and the relative phase between them was 
initially undetermined. 

One way to initialize a system with a well defined 
phase is to construct the analogue of a beam splitter 
whereby a single condensate is dynamically split into 
two coherent parts, which serve as the two "arms" of 
the interfcrometeiiSi^. In recent experiments this was 
achieved by raising a potential barrier along the axis of a 
quasi one-dimensional condensate^. The split is applied 
slowly compared to the transverse frequency of the trap, 
but fast compared to the longitudinal time scales. Thus 
each atom is transferred in this process to the symmetric 
superposition between the two traps without significantly 
changing its longitudinal position. An illustration of the 
procedure is given in Fig. ^ In repeated experiments 
the condensates are released from the trap to probe the 



interference at varying times after the split. Immedi- 
ately following the split the two condensates are almost 
perfectly in phase. However repeated measurements at 
longer times show that the phase distribution becomes 
gradually broader until it becomes uniformly distributed 
on the interval [0, 27r]. The accuracy of interferometric 
measurements is limited by such dephasing. 

In this paper we develop a theory of dephasing in in- 
terferometers made of one dimensional Bose gases. We 
assume that such systems can be well isolated from the 
noisy environment, so that this source of decoherencsiS 
is set at bay. Another dephasing mechanism that has 
been discussed extensively in the context of condensates 
in double wells is "quantum phase difFusion"iii^. The 
uncertainty in the particle numbers in each well brought 
about by the split, entails a concomitant uncertainty in 
the chemical potential difference, which leads to broad- 
ening of the global relative phase on a time scale of 
T£) ^ h^/N / ^. Here N is the total particle number and /i 
is the chemical potential. Note that this is essentially a fi- 
nite size effect, td diverges in the proper thermodynamic 
limit in which the trap frequency is taken to zero while 
the density is kept constant. Here we show that dephas- 
ing in quasi one dimensional systems is dominated by 
a new mechanism that involves quantum fiuctuations of 
the local relative phase field rather than the global phase 
difference between the condensates. In contrast to the 
usual phase diffusion this is a bulk mechanism that sur- 
vives the thermodynamic limit. Indeed, it is well known 
that in one dimensional liquids quantum phase fiuctua- 
tions arc extremely effective and prevent the formation 
of long range order. Here we show how these fluctua- 
tions destroy long range order in the relative phase, that 
was imposed on the system by the initial conditions. We 
account for the phase fiuctuations using a quantum hy- 
drodynamic description of the Bose liquidsii. The phase 
coherence between the condensates is then shown to de- 
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FIG. 1: Interferometer setup. At time t — a, single one-dimensional condensate is split into two coherent parts by raising 
a high potential barrier along the condensate axis. After a wait time t the condensates are released from the trap and let to 
expand. After a sufficient expansion time the cloud is imaged by a probe beam sent along the condensate axis and the phase 
of the resulting interference pattern is recorded. This process is repeated many times in order to obtain a phase distribution 
at a number of different wait times t. 



cay exponentially in time. The dephasing time we obtain 
is a simple function of intrinsic parameters: the dimcn- 
sionlcss interaction strength, the density, and the sound 
velocity. These predictions are verified and extended to 
finite temperatures, using direct simulation of the dy- 
namics in the microscopic hamiltonian. The results are 
also compared to the dephasing time measured in Ref. 
0, and are found to be in quantitative agreement. 

It is interesting to point out the essential difference 
between the problem we consider and dephasing of sin- 
gle particle interference effects as considered for exam- 
ple in mesoscopic electron systems. As in our problem, 
phase coherence in mesoscopic systems can be defined 
by an interferometric measuremen^i^. In a Fabry-Perot 
setup, oscillations of the conductivity as a function of 
applied Aharonov-Bohm flux vanish exponentially with 
system size, thereby defining a characteristic dephasing 
length. Because what is measured is ultimately the DC 
conductivity, the problem may be recast in terms of lin- 
ear response theorjsi^. By contrast, dynamic splitting 
of the condensate in the ultra-cold atom interferometer 
takes the system far from equilibrium and the question 
of phase coherence is then essentially one of quantum dy- 
namics. The system is prepared in an initial state deter- 
mined by the ground state of a single condensate, which 
then evolves under the influence of a completely different 
hamiltonian. that of the split system. Dephasing, from 
this point of view, is the process that takes the system 
to a new steady (or quasi-steady) state. 



II. HYDRODYNAMIC THEORY 

We start our analysis by considering the hydrodynamic 
theory that describes low energy properties of one dimen- 
sional Bose liquidsii. The hydrodynamic hamiltonian for 
a pair of decoupled condensates is that of two decoupled 



Luttinger liquids (we set h = 1 throughout) 



TT K 



(1) 



where c is the sound velocity and K is the Luttinger pa- 
rameter that determines the decay of correlations at long 
distance. IIq, is the density fluctuation operator conju- 
gate to the phase ipa ■ The smooth component of the Bose 
field operator is given by ipa ~ y^'^^Pi^ 'Pa), where p is 
the average densitjiii. 

The operator corresponding to the interference signal 
between the two condensates is given byi^: 



A: 



P 



dx 



(2) 



where L is the imaging length. For a pair of decoupled 
condensates at equilibrium, {A) = while {A^) > Oi^. 
Therefore the interference pattern displayed in repeated 
experiments has a finite amplitude but its phase is com- 
pletely random. On the contrary, in this work we consider 
a pair of condensates that are prepared out of equilib- 
rium, with a well defined relative phase. In this case (A) 
is expected to be non vanishing at the time of the split 
and decay in time as the fluctuation in relative phase 
grows. 

Calculation of the time evolution of (A) is greatly 
simplifled by the fact that the hydrodynamic theory is 
quadratic. In particular, this allows decoupling in the 
hamiltonian of " center of mass" and relative phase flelds 
(p+ = {ipi + ip2)/'2 and (fi- = Lpi~ if2- In other words, the 
hamiltonian can be rewritten as a sum of two commuting 
harmonic terms H = i7+((/3+,n+) + iJ_ ((^_ , n_), where 
n± are the conjugate momenta of ip±. Moreover, the 
splitting process described above ensures that the ini- 
tial state may be factorized as ^[(/J-i-, = 'I'+[((9+] x 
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_[(/?_]. Because during the split all atoms are sim- 
ply transferred to the symmetric superposition without 
changing their axial state, the wave function is 
identical to the wave function of the single condensate 
before the split. If this condensate is at finite tempera- 
ture \E'+ is replaced by the appropriate density matrix. 
On the other hand is determined by the splitting 
process and will generally be strongly localized around 
f-{x) = 0. In fact we can find the approximate form of 
this wavefunction, which will be an important input for 
the time dependant calculation. 

Consider a region of size ^/j, which is the length scale 
on which the hydrodynamic variables are defined. Let 
ni and n2 be the operators corresponding to the par- 
ticle number in each of the condensates within this re- 
gion. The splitting process delocalizes the particles be- 
tween the two condensates, resulting in a roughly gaus- 
sian distribution of the relative particle number ?i_ = 
(ni— n2)/2, with an uncertainty of order (rii + 71,2) /4 = 

pS,h/'2- The phase ip- defined on this grain in the 
hydrodynamic theory is canonically conjugate to n_. 
Therefore the wave function in the phase representa- 
tion is approximately gaussian within {ip{x)ip{x')) ~ 
S{x-x')/{2p). That is 



\Vq\ 



(3) 



From here on we drop the subscript (-) where no ambi- 
guity may arise. 

We are now ready to compute (A) which depends only 
on if- . The wave-function [ip-] evolves in time under 
the influence of the harmonic hamiltonian 



2n , Kq^ , 



The wavefunction remains gaussian at all times: 
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where r]q(ipq, ip*, t) is a pure phase and 

Wq{t) = 



1 ^ 2tt'^ ps\T?(cqt) 



(6) 



The expectation value {A) = ^ in the gaussian wave 
function jSJl is given by 



A(t) = p(e"^(^)) = pe-^^o ^l-^'l'^ 



pe 



-Bit) 



(7) 



In the thermodynamic limit [L ^ 00) the summation 
in lO can be converted to an integral to obtain: 



9{t) 



2-K 



90 



dq 



1 9, , 2'!T^psm'^(cqt) 
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(8) 



The last equality is valid at times t > ^/i/c, when the inte- 
grals become independent of the high momentum cutoff. 
In addition we require that the particle number in a heal- 
ing length S,h is large, i.e. p^h 3> 1, this is always satisfied 
at weak coupling K ^ 1. Thus we find that the phase 
coherence decays exponentially: 



-t/7 



Ait) = Aoe 



with the characteristic time 



T = — {K/ttY 

cp 



(9) 



(10) 



Note that the dephasing time t diverges, as it must, in 
the non interacting limit (K — > 00 and c — *■ 0). Clearly 
the initial state induced by the split is an eigenstate of 
non interacting particles and therefore all observables in- 
cluding the coherence must be time independent in this 
limit. 

For a finite system, the discreteness of the sum in O 
must be taken into account at times t ^ L/c. In such 
a case, the dominant contribution to g{t) arises from the 
q = term. This term yields a ballistic broadening of 
the phase ((/Jq^g) oc i^, which results in a gaussian decay 
of the coherence: A{t) = ^oexp(— t^/r|) where tl = 
v/rL/2c. 

It is interesting to note that one can also express the 
time scale associated with the finite size as oc a/ZV /p. 
Written in this way, it is clear that the gaussian term 
is just the "quantum phase diffusion" discussed in the 
context of double well systems^. It stems from the un- 
certainty relation between the uniform component of the 
relative phase and the difference in total particle num- 
bers between the wells. The uniform component of the 
interaction, which drives the ballistic broadening of ipq=o 
scales like l/L, making the gaussian term irrelevant in 
large systems. By contrast, the exponential dephasing 
due to the internal modes is length independent. The fact 
that these modes produce exponential dephasing is spe- 
cial to one dimensional interferometers. A similar analy- 
sis of two parallel planar condensates yields logarithmic 
divergence in ((SJ, which implies power law decay of the 
coherence as A oc {(^h/ct)". We note that in this case K 
is not dimensionless. The power a is non universal and 
depends on the short distance cutoff. In three dimensions 
on the other hand the only contribution to dephasing is 
the usual phase diffusion due to finite size and there is 
no bulk contribution. 

It is striking that the hydrodynamic theory developed 
in this section predicts a dephasing time which is inde- 
pendent of temperature. Temperature enters the initial 
condition only through the density matrix for the sym- 
metric degrees of freedom, whereas the phase coherence 
between the condensates depends on the evolution of the 
relative phase. Because the symmetric and antisymmet- 
ric fields are completely decoupled within the harmonic 
theory, the temperature associated with the symmetric 
degrees of freedom does not affect the dephasing. How- 
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ever, it is clear that the fuh microscopic hamiltonian con- 
tains anharmonic terms that do couple those degrees of 
freedom. At equilibrium the non linear terms are irrel- 
evant (in the renormalization group sense) and do not 
affect the asymptotic long wavelength low energy corre- 
lations. However it is natural to question the validity 
of the hydrodynamic description for computing time de- 
pendent properties. In particular, here the system is pre- 
pared out of equilibrium in a state with extensive energy 
relative to the ground state. Then it is a priori unclear 
that the low energy correlations given by the hydrody- 
namic theory are sufficient to describe the dynamics. In 
the next section we shall test the predictions of the hy- 
drodynamic theory and obtain temperature dependent 
corrections to the dephasing by computing the dynamics 
within a microscopic model of the twin condensates. 

III. NUMERICAL RESULTS 

For the purpose of numerical calculations we shall con- 
sider a lattice hamiltonian, the Bose-Hubbard model on 
twin chains 

2 jj 
H = ^(^LWl + H.c.) + Y ^ ^mK. - 1). 

i a — 1 Qz 

(11) 

Here creates a Boson on site i of chain a. The model 
(111(1 can also describe continuum systems, such as the 
one in Ref. 0, if the average site occupation n = {rij) is 
much less than unity. 

As mentioned above, we are interested in the time evo- 
lution of the expectation value A{t) = {b\^{t)b2i{t)) . 
The average is taken over the wave function or density 
matrix at the time of the split. We will calculate A{t) us- 
ing a semiclassical approach called the truncated Wigner 
approximation (TWA, cf. Ref. Let us first briefly 

review the strictly classical approximation to the dynam- 
ics. The usual procedure is to write the Heisenberg equa- 
tions of motion for hi and b\ and then replace these op- 
erators by complex classical fields ■0^ and if)*. This leads 
to the lattice Gross-Pitaevskii (GP) equations 

= -J{i^a,i~l + 1pa.i+l) + ?7|V'a,jP'0a,i- (12) 

at 

Given the initial condition '0Q,i(O) one can integrate the 
GP equations to find the value of the fields at any time t 
and obtain Aci{t) = X)i ^ii(*)''/'2i(t)- If the split is fully 
coherent ■(/'ii(O) = V'2i(0). Since the evolution in the two 
chains is described by identical equations, this equality 
persists to all subsequent times. Thus in the absence of 
external noise sources, the classical dynamics cannot ac- 
count for dephasing. By contrast, quantum fluctuations 
provide an intrinsic dephasing mechanism. 

Quantum corrections modify the dynamics in two 
ways. First, they introduce "quantum noise" to the ini- 
tial conditions. The fields tpai and ■0*^ originate from non 



commuting quantum operators bai and b\^ which cannot 
be determined simultaneously. Therefore the unique clas- 
sical initial condition should be replaced by an ensemble 
of initial conditions characterized by a quantum distri- 
bution. Second, the classical trajectories determined by 
(|12|1 are supplemented by additional quantum paths. 

It can be shown that the leading quantum correction 
to the dynamics enters through the initial conditions'^. 
The essence of TWA is to integrate the GP equations 
starting from an ensemble of initial conditions, which are 
given quantum weights derived from the initial density 
matrix. Thus within this approximation: 

A{t)^ Jv{rja,i'jo.}opMjja,Mo)Acl{t). (13) 

Here 

drfdrj {i^ ~ -\ Po\ + ^) 
^^-M'-iM\Uv'^-vr)^ (14) 

is the Wigner representation of the initial density matrix 
of the system and | tp ) denotes a coherent state with 
eigenvalue t/j of the boson annihilation operator. We note 
that pw should not be thought of as a probability distri- 
bution, because in general it can assume negative values. 

The main difficulty in applying the recipe H13|) to com- 
pute A{t) is the need to find the Wigner distribution of 
the split condensate at the time of the split. To overcome 
this problem we use the following procedure. Rather than 
tackling a split interacting system we start the calcula- 
tion from a single non-interacting condensate, where we 
have an exact expression for pviT^- The evolution of the 
Wigner distribution from a non-interacting system to an 
interacting one is done using the TWA by slowly increas- 
ing the interaction constant, U, from zero to the desired 
value (see footnote 0|). The heating induced by the 
time dependant hamiltonian is controlled by the rate at 
which U is increased. The next step is the split of the 
condensate. It involves the doubling of the degrees of 
freedom at each site. Because of the way the split is car- 
ried out, the field ipi of the single condensate is simply 
copied to the symmetric field ip+^i = (■0ii-|--02i)/\/2. The 
new degree of freedom ip-,i is chosen at this point from 
a Wigner distribution 

p^-\{r-,^-}o) = n (l) e-'I'^-'l' (15) 

which corresponds to the vacuum of the boson b^j ~ 
{bii — bii)l\pl. The calculation is continued using the 
dynamics of the twin condensates. 

The condensate (prior to the split) is prepared at var- 
ious temperatures using two different methods: (i) By 
initializing the non interacting condensate in a finite tem- 
perature density matrix, (ii) By varying the rate by 
which the interactions are switched on to induce a con- 
trolled amount of heating. The temperature of the con- 
densate before the split is parameterized by measuring 
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FIG. 2: Decay of the coherence calculated with the truncated 
Wigner approximation for U = 1, J = 1 and n = 7.5 com- 
pared to the result of the hydrodynamic theory (dots). The 
time axis is scaled to the hydrodynamic dephasing time tq 
llOI I. Lines correspond to different initial temperatures pa- 
rameterized by the correlation length that existed before 
the split. 
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and those of the hydrodynamic theory l^^^- pBi- 
ticular the value of the Luttinger parameter is given by 
K w -!Ty/2Jn/U , and the sound velocity is c = \/2JnU . 
The short distance cutoff of the hydrodynamic theory 
may also be extracted from 1)11(1 . by finding the length 
over which amphtude fluctuations decay. It is given byiS 
Ih = \/26i ~ ayjAjJUfL. Increasing n decreases this 
length scale until it reaches its lower bound which is 
one lattice constant a. When Ih ^ a the model Hll(l 
in effect describes a continuum system, while for a large 
site occupancy ^ a it describes a lattice of coupled 
condensates. We perform simulations in both regimes. 
The continuum regime is much more demanding com- 
putationally in the sense that attaining a given thermal 
correlation length in this regime requires a much slower 
activation rate of the interaction. This limited the tem- 
perature range we could study in the continuum regime. 

The numerical results are summarized in Fig. [3 and 
Fig. 13 The dephasing time, r, is extracted by fitting 
A{t) to Eq. 0. At low temperatures, t is in good agree- 
ment with the hydrodynamic result (|10() and shows very 
weak temperature dependence. On the other hand above 
a crossover temperature T*, the dephasing time gains 
significant temperature dependence. The crossover scale 
is set by the condition ^t* ~ Ih- This is not surprising, 
given that the hydrodynamic theory is expected to break 
down when correlations decay over a length scale shorter 
than its short distance cutoff. The agreement between 
the numerical results and the hydrodynamic theory at 
low temperatures is highly non-trivial in view of the fact 
that the two approaches are based on entirely different 
approximations. Most importantly, the temperature in- 
dependent result H10() relied on the decoupling of sym- 
metric and antisymmetric degrees of freedom within the 
harmonic theory (Q. Such decoupling does not exist in 
the microscopic model (|llll nor in the TWA dynamics. 
In the next section we shall compare the results of the 
hydrodynamic theory to experimental data. 



IV. COMPARISON WITH EXPERIMENTS 



FIG. 3: Dephasing time r as a function of the thermal corre- 
lation length, which is a measure of the inverse temperature 
before the split. Data is taken in different regimes of micro- 
scopic parameters. All data is seen to collapse on approxi- 
mately the same curve when the dephasing time is scaled to 
the theoretical time scale and the thermal correlation 

length scaled by Ih (see text). 



the correlation length oc 1/T associated with the de- 
cay of the correlation function (fej&i+r)- We verified that 
the final result for A{t) depends only on and not on 
the method used to achieve finite temperature. 

Finally, in order to compare the numerical results to 
the prediction of the hydrodynamic theory we use the 
fact that at weak coupling there are simple relations 
between the parameters of the microscopic model (|11|) 



The dephasing time ((10|1 is written in terms of param- 
eters of the hydrodynamic theory. In order to compare 
with experiments we have to translate these parameters 
to numbers that are relevant to a specific experimental 
situation. For this purpose we will consider the setup of 
the experiments described in Ref. 0. 

This system consists of bosons with contact inter- 
actions parameterized by a dimensionless interaction 
strength 7. At weak coupling, 7^1, the following rela- 
tions holcii^: K ~ 7i'/-y/7 and c = h{p/M)y/j, where M 
is the mass of an atom. Substituting these into (|10|l we 
obtain: 



/2M 



7 



3/2 ■ 



(16) 



Now consider a one dimensional tube geometry with 
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transverse trap frequency U!±. If the oscihator length 
lj_ is much larger than the s — wave scattering length a^, 
as is the case in Ref . , we can use the approximate ex- 
pression 7 w 2Muj±as/Hp^^^. This allows us to write 
the dephasing time using parameters easily determined 
in the experiment: 

" = (^) ^''^ 

Note that p is the density per condensate in the split 
system, that is half the density of the initial condensate. 

Taking the scattering length of rubidium-87 atoms 
flg = 105ao, the transverse trap frequency lu± ~ 2tt x 
2.1 kHz and density of about 50 atoms per micronS^, 
we obtain a dephasing time r w 4.3ms. To estimate the 
dephasing time in the experiment^ we use the data for 
the phase broadening assuming a von Mises distribution 
of the phase {f{9) = exp(KCOS0)/27r/o(K)). This gives 
approximately 2ms, slightly shorter than our theoretical 
estimate. We note that data points marked with squares 
in Fig. ^correspond to parameters relevant to the experi- 
ment. Unfortunately the temperature in that experiment 
is not known to a good accuracy. In the future it would 
be interesting to look for the universal temperature de- 
pendence implied by the numerical results. 

V. DISCUSSION AND CONCLUSIONS 

It is well known that quantum fluctuations prevent 
long range phase order from forming in one-dimensional 
Bose liquids. This phenomenon is most conveniently de- 
scribed within the framework of the Luttinger liquid or 
hydrodynamic theory. Here we used this framework in a 
non equilibrium situation to show how quantum fluctua- 
tions destroy long range order that was imposed on the 
system as an initial condition. The outcome is a simple 
formula describing exponential decay of phase coherence 
in interferometers made of one-dimensional condensates. 
The dephasing time found in this way provides a fun- 
damental limit on the accuracy of such interferometers. 
We did not discuss the situation in which the interfer- 
ometer is prepared in a number squeezed initial state. 
Such an initial condition would slow the usual phase dif- 
fusion procesai^ and is expected to similarly affect the 
bulk mechanism discussed in this paper. 

Interestingly the dephasing time measured in Ref. 
is only slightly shorter than the prediction of the hydro- 
dynamic theory. We therefore conclude that quantum 
phase fluctuations were probably a dominant dephasing 
mechanism in that experiment. 



The validity of the Luttinger liquid framework out of 
equilibrium is not ensured a priori. We therefore test its 
predictions against numerical simulations of the micro- 
scopic model using the Truncated Wigner approximation. 
At low temperatures, the simulation results were essen- 
tially temperature independent and in good agreement 
with the hydrodynamic theory. On the other hand at 
temperatures above a crossover scale set by the healing 
length, the dephasing time displayed considerable tem- 
perature dependence. Indeed the Luttinger liquid theory 
is expected to break down above this temperature scale. 

Phase fluctuations are expected to be weaker at higher 
dimensions. For example three dimensional condensates 
can support long range order even at flnite temperatures. 
In this case the long range order in the relative phase im- 
posed as an initial condition of the interferometer is resis- 
tant to the phase fluctuations. The situation with planar 
condensates is more delicate, two dimensional systems 
can in principle sustain long range order at equilibrium 
at zero temperature. However the initial condition of the 
interferometer drives the system out of equilibrium. The 
hydrodynamic theory yields power law dephasing in this 
case. However the power is non universal and further 
work is needed to test the validity of the hydrodynamic 
theory. 

Finally we would like to point out that the dephas- 
ing process considered in this paper is a mechanism that 
brings the system, from a non equilibrium state imposed 
by the initial conditions to a new steady state. According 
to the hydrodynamic theory this steady (or quasi-steady) 
state is not yet thermal equilibrium. In particular we find 
that the off diagonal correlations along the condensates 
are distinctly non thermal at steady state^i. It is an 
interesting question whether there are processes, much 
slower than dephasing, that eventually take the system 
toward thermal equilibrium. This issue can be addressed 
by experiments. After the phase has randomized, cor- 
relations along the condensates can be measured using 
analysis methods developed in Refs. [323 1. Thus, we 
propose that interferometric experiments can serve as de- 
tailed probes to address fundamental questions in non 
equilibrium quantum dynamics, supplementing measure- 
ments of global properties previously used to touch on 
these issues^. 
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